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Abstract 

We propose a simple explanation for the connection between chiral 
symmetry restoration and deconfinement in QCD at high tempera- 
ture. In the Higgs description of the QCD vacuum both spontaneous 
chiral symmetry breaking and efi'ective gluon masses are generated by 
the condensate of a color octet quark-antiquark pair. The transition 
to the high temperature state proceeds by the melting of this con- 
densate. Quarks and gluons become (approximately) massless at the 
same critical temperature. For instanton-dominated effective multi- 
quark interactions and three light quarks with equal mass we find a 
first order phase transition at a critical temperature around 170 MeV. 
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1 Octet melting 



Lattice simulations for QCD at high temperature [|I| find empirically that 
chiral symmetry restoration and deconfinement occur at the same critical 
temperature Tc Sufficiently below Tc the equation of state is rather well 
approximated by a gas of mesons. Gluons and quarks or baryons play no 
important role. Above the dominant thermodynamic degrees of freedom 
are gluons and quarks. The change in the relevant number of degrees of 
freedom occurs rather rapidly and may be associated with a thermodynamic 
phase transition. 

The transition to the high-temperature state of QCD is often called "de- 
confinement transition" since the gluon modes with momenta ~ (^^)^ 
characteristic for a thermal state behave at high temperature close to a gas 
of free particles. (For experimentally accessible temperatures above the phase 
transition substantial corrections to the equation of state of a relativistic gas 
remain present.) This also holds for the fermionic degrees of freedom which 
appear essentially as free quarks. Having a closer look the picture of weakly 
interacting gluons is actually somewhat oversimplified. In fact, the long- 
distance behavior of the time-averaged correlation functions in high temper- 
ature QCD corresponds to a "confining" three-dimensional effective theory 
[^] with a strong effective coupling and a temperature-dependent "confine- 
ment scale" [0. Furthermore, the transition between confinement and 
weakly interacting "free" quarks is less dramatic than it may seem at first 
sight. In particular, it needs not to be associated with a true phase tran- 
sition. In presence of light quarks, the temperature-induced changes in the 
long-distance behavior of the heavy quark potential are only quantitative [0]. 
String breaking occurs for arbitrary temperature. 

Nevertheless, one expects a temperature range where the properties of 
strong interactions change rapidly, refiecting "deconfinement" . We associate 
here "deconfinement" with the transition to a thermal equilibrium state for 
which (1) the bosonic contribution to the free energy can be approximated 
by the one of a free gas of eight massless spin-one bosons and, correspond- 
ingly, (2) possible masslike terms in the effective "gluon-propagator" are of 
the order of the temperature or smaller. We do not deal in this paper with 
temperature-dependent changes in the heavy quark potential and tempera- 
ture effects in the dynamics of string breaking. 

On the other hand, the chiral properties of QCD are also expected to 
undergo a rapid change at some temperature. The vacuum is characterized 
by a chiral condensate of quark-antiquark pairs. This order is destroyed 
at high temperature and chiral symmetry is restored. For vanishing quark 
masses this implies a true phase transition. In this case chiral symmetry 
restoration at high temperature is signalled by a vanishing order parameter. 
At high temperature the massless quarks contribute to the free energy as 
a relativistic gas with 3A^/ fermions whereas below a critical temperature 
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the fermionic degrees of freedom appear as massive baryons without much 
relevance for thermodynamics. This quahtative change in the fermionic con- 
tribution remains true also in presence of small nonvanishing quark masses. 
We associate the chiral transition at high T with a rapid decrease of jump 
in the chiral order parameter and the appearance of a relativistic quark gas 
contribution to the free energy. 

The qualitative feature that at high temperature deconfinement and chi- 
ral symmetry restoration come in pair is not much of a mystery: modes 
with momenta ~ T dominate the thermal state and the interactions be- 
tween these modes become comparatively weak. Symmetry restoration at 
high T is a common phenomenon. It is a longstanding puzzle, however, 
why deconfinement and chiral symmetry restoration occur quantitatively at 
the same temperature. Analytical approaches that have been investigated 
so far either concentrate on the quark degrees of freedom - for example in 
the context of Nambu-Jona-Lasinio models This reflects important as- 
pects of the chiral transition but fails to describe the effects of confinement 
and the thermodynamic equation of state of QCD. Other models deal with 
pure QCD without including the important aspects connected with the pres- 
ence of light quarks. So far we are not aware of any successful quantitative 
analytical description accounting simultaneously for the deconfinement and 
chiral symmetry restoration aspects of the QCD phase transition. In view 
of the vaste experimental efforts and the important progress in numerical 
simulations [|10| , even an oversimplified analytical description of the most 
salient characteristics of this transition would be very welcome. 

One of the main difficulties for any analytical description are the differ- 
ent relevant degrees of freedom above and below the critical temperature. 
Whereas for high T a gas of quarks and gluons becomes a reasonable ap- 
proximation, the low temperature physics is described by an interacting pion 
gas with far less degrees of freedom. The fermionic degrees of freedom at low 
T are baryons which are Boltzmann-supressed. Any quantitative analytical 
description must be able to describe both gluons and pions as well as quarks 
and baryons simultaneously. This is mandatory at least in the vicinity of the 
critical temperature which is precisely characterized by the transition from 
one set of dominant degrees of freedom to another. The recently proposed 



Higgs-picture of the QCD-vacuum |TT|, |1^ offers such a possibility since all 



relevant degrees of freedom are described at once. A reasonable picture of 
the QCD-phase transition may also serve as a test for these ideas. 

In this paper we propose a simple explanation of the simultaneous occur- 
rence of the deconfinement and chiral transition and present first quantitative 
estimates. Our explanation is based on the proposed new understanding of 
the QCD-vacuum by "spontaneous breaking of color" or gluon-meson du- 
ality |rT|. In this picture of the vacuum the condensate < % > of a color 
octet quark-antiquark pair leads to spontaneous breaking of the color sym- 
metry. (As for the electroweak standard model there exists an equivalent 



2 



gauge invariant formulation.) The gluons acquire a mass ~< x > by the 
Higgs mechanism and become integer charged. Gluon-meson duahty asso- 
ciates the massive gluons with the physical vector mesons p, K* and uo. This 
is possible since they carry the appropriate integer electric charges as well 
as isospin and strangeness. Also the quarks carry integer charges after spon- 
taneous color symmetry breaking and can be associated with the baryons 
(p, n, A, S, H) plus a heavy singlet. They become massive due to the sponta- 
neous breaking of the chiral flavor symmetry by < x > and a similar singlet 
qq- condensate. This association between quark fields and baryons is called 
quark-baryon duality. It is the analogue of color-flavor locking |]l3l for a high 
baryon density. 

It has been argued on phenomenological grounds |]ll| that the dominant 
contribution to spontaneous chiral symmetry breaking arises from the octet 
condensate < x >. The same condensate therefore explains both confinement 
and spontaneous chiral symmetry breaking. In this description the solution to 
the puzzle why deconfinement and chiral symmetry restoration are connected 
becomes obvious. At some critical temperature Tc the value of the octet 
condensate < x> drops rapidly. At this temperature the mass of the gluons 
therefore drastically decreases and the deconfinement transition happens. 
In the limit of three massless quarks the expectation value < x > exactly 
vanishes at T^, reflecting chiral symmetry restoration. Typically, the color- 
octet and -singlet gg-expectation values influence each other. In particular, 
an octet condensate always induces a singlet condensate. It seems plausible 
that both the octet and singlet condensates vanish simultaneously for T > T^. 
Then chiral symmetry gets completely restored for T > T^, leading to a true 
phase transition characterized by the order parameter of chiral symmetry 
breaking. The critical temperature for the phase transition is the "melting 
temperature" for the octet condensate. The "deconfinement temperature" 
and "chiral symmetry restoration temperature" are therefore identical. As a 
consequence of the change in the order parameter the quantum numbers of 
the excitations in the "hadronic phase" for T < Tc differ from the ones in 
the "quark-gluon phase" for T > Tc. 

Quark mass effects modify the details of the transition. For a more quan- 
titative description of "octet melting" we omit in this paper the difference 
between the current quark masses and consider three light quark flavors 
with equal mass. For our quantitative analysis we will need information 
about the condensates in the vacuum which are related to the form of the 
effective quark interactions. We investigate a rather wide class of effective 
multiquark interactions induced by instanton effects. In the chiral limit of 
three massless quarks we find a first-order transition with critical tempera- 
ture Tc ~ 130 — 160 MeV. In this limit the mass of the eight pseudoscalar 
Goldstone bosons {TT,K,ri) vanishes at low temperature, mps = 0. For real- 
istic current quark masses corresponding to a nonzero average pseudoscalar 
ip5 ~ i^M]^ + M^)/3, we again get a first-order transition with a 
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somewhat larger Tc ~ 170-180 MeV. 

In this note we do not deal explicitly with the color-singlet quark-antiquark 
pair which plays a subdominant role. The relevant thermodynamic potential 
is then given by the temperature-dependent effective potential for the octet 
condensate 

[/(x,T) = f/o(x) + Af/(x,T) (1) 

Here Uq encodes in a bosonic language for scalar gg-composites the informa- 
tion about the multiquark interactions in the vacuum, whereas Af/ accounts 
for the thermal fluctuations. We use a scalar field Xij,ab ~ '4'jbi^ai — l'4'kb4'akSij 
for the octet quark-antiquark pair with i,j = 1...3 color indices and a,b = 
1---Nf flavor indices for massless (or light) quarks. For Nf = 3 it is sufficient 
to evaluate the effective potential for the direction of the condensate [|Tl|] 

1 1 

< Xij,ab >= -^Xi^iaSjb - -SijSab) (2) 

In the limit of three equal quark masses this condensates preserves a vector- 
like S'f/(3)-symmetry which contains the generators of isospin and strangeness. 

We emphasize that in our picture the thermodynamic quantities depend 
on the number of light flavors in an important way. For example, the vacuum 
condensates for two-flavor QCD {Nf = 2) are discussed in ||l^ and differ 



substantially from the three-flavor case. Needless to say that in absence of 
light quarks the chiral aspects of the phase transition are completely different. 
We restrict our discussion here to the realistic case Nf = 3. Quark mass 
effects are taken into account except for the S't/(3)-violation due to the mass 
differences. The thermodynamic quantities of interest can then be extracted 
from the behavior of the minima of U{x, T). For T = the minimum occurs 
for a nonvanishing octet condensate Xo 7^ 0- As the temperature increases, 
the thermal fluctuations induce a new local minimum at the origin x = 0- At 
the critical temperature the two minima are degenerate. The discontinuity 
characteristic for a flrst-order phase transition corresponds to the jump to 
the absolute minimum at x = for T > Tc. 

The remainder of this paper is devoted to a more detailed estimate of 
U{x,T). We will employ a mean fleld-type calculation which takes into ac- 
count, nevertheless, the most important effects of the strong interactions. We 
will devote some care to the proper implementation of the mean fleld calcu- 
lation despite the fact that some other approximations in our approach (like 
the effective vacuum potential) remain rather crude. The reason is that we 
want to describe correctly the interacting pion gas at low temperature within 
a linear formulation. The linear formulation is necessary for a description 
of symmetry restoration at Tc, since for T > the Goldstone bosons are 
combined with the a -scalar into linear multiplets of the flavor group. On 
the other hand, the nonlinear Goldstone excitations dominate the free energy 
at low T and chiral perturbation theory gives a valid description. We will 
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see that a linear mean field description easily leads to incorrect results for 
the nonlinear excitations unless done with sufficient care. Furthermore, we 
do not want to obscure the possibility of a second order phase transition by 
using a method which is too rough. Again, unless the bosonic degrees of 
freedom are treated with care, a mean field computation can easily produce 
a spurious first order jump even in a situation of a second order transition 
or a crossover. This is well known from the study of scalar field theories or 
the electroweak phase transition. 

Our paper is organized as follows: Sect. 2 addresses the form of the 
vacuum potential Uo{x), which was estimated previously fl^ by an instanton 



calculation. In sect. 3 we discuss the dependence of the masses of the most 
relevant excitations on the octet condensate x- This is the basis for our 
mean field-type calculation of the temperature effects AU{x,T) in sect. 4. 
A first summary of quantitative results for the chiral phase transition is 
given in sect. 5 in the chiral limit of massless quarks. For the instanton- 
induced potential Uo{x) proposed in sect. 2 we find a critical temperature 
Tc=154 MeV, consistent with lattice estimates [0]. Sect. 6 is devoted to 
the high-temperature phase and the equation of state. Sufficiently above Tc 
one observes a weakly interacting gas of gluons and quarks. In sect. 7 we 
present analytical estimates for the chiral phase transition. In particular, for 
three massless quarks the critical temperature is related to the mass of the 
?7'-meson and its decay constant fe by 

= O.OISMJ (3) 
with i? < 1 a constant of order one. 

We turn to the value of the chiral symmetry-breaking order parameter in 
the low temperature phase in sect. 8. This needs a meaningful treatment of 
the Goldstone boson fluctuations which we propose in sect. 9. In sect. 10 we 
verify that our prescriptions are consistent with chiral perturbation theory 
for low enough T. Sect. 11 discusses the effects of nonvanishing current 
quark masses. We consider equal current quark masses which correspond 
to an average mass of the pseudoscalar Goldstone bosons Mps =390 MeV. 
Up to S'?7(3)-violating quark mass differences this corresponds to a realistic 
situation. We find Tc =170 MeV, again consistent with lattice simulations 
[0]. In the low temperature phase the Higgs contribution to the mass of the 
vector mesons decreases substantially as the temperature increases. In the 
vicinity of we find a Higgs contribution to the average vector meson mass 
^ 300 MeV and to the average mass for the light baryons ~ 600 MeV. Those 
have to be supplemented by thermal mass contributions. A decrease of the 
vector meson mass [|T5| may be relevant for the dilepton spectrum observed 



in heavy ion collisions. Finally, we investigate in sect. 12 a possible lattice 
test of our scenario by establishing the existence of stable Za-vortices in the 
low temperature phase. Sect. 13 summarizes our conclusions. 
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2 Vacuum-effective potential 



The effect of the thermal fluctuations has to be compared with the vacuum- 
effective potential t/o(x) for the octet. We will not need here many details of 
the precise shape of Uq. The most important parameter for the determination 
of the critical temperature Tc for the phase transition will turn out to be 
the difference Uo{0) — f/o(Xo)) with xo the vacuum expectation value of the 
octet. To be specific, we concentrate here on an instanton-induced potential 
which is motivated by the observation that instanton effects are a plausible 
candidate for the dynamics that lead to a non- vanishing octet condensate [|1^ . 
Beyond the octet gg-condensate realistic QCD exhibits also a nonvanishing 
color-singlet gg-condensate a. In order to concentrate our discussion on the 
essential aspects, we assume here that a is "integrated out" by solving its 
effective field equation^] as a functional of x- We take here a ~ x (i^i the 
relevant range) such that all qq condensates are accounted for by x- 

For Nf = 3 and small values of qq the instanton effects lead to a cubic 
term, Uan ~ X'^ ~ (qq)^- For large values of x the nonvanishing gluon mass ~ 
X leads to an effective infrared cutoff for the instanton interactions, resulting 
in Uan ~ X~^^ 101 • A first quantitative estimate of the instanton-induced 



effective potential can be found in |jTj] and we exploit its consequences in 



the present work. However, in view of the uncertainties of this estimate we 
will be satisfied here with a simplified version which respects the asymptotic 
properties for small and large x 

(A) : Uoix) = A{ - ^i?.„Xox'(l + ^i-YT' + 

11 11 Xo 

(l-i?.„)(x'-2x^x')+X^} (4) 

The instanton contribution is the one ~ Ran and we have added another 
contribution ~ (1 — Ran) which arises from non-anomalous interactions not 
related to instantons, as, for example, multiquark interactions mediated by 
gluon exchange. Presumably the vacuum potential is dominated by instanton 
effects such that Ran is expected close to one. In our conventions x is 
real and we have parametrized the potential such that the absolute minimum 
of Uo occurs at xo, with f/o(x = 0) = Xxq, A > 0. The additive constant is 
fixed such that the pressure vanishes in the vacuum. In absence of quark 
masses this implies Uo{xo) = 0. 

In order to investigate how sensitively our results depend on the precise 
shape of the potential, we consider for Nf = 3 also a potential which only 
incorporates the breaking of the discrete symmetry x ~^ —X by the anomaly, 



^The effective kinetic term and interactions of the ^-field include therefore contributions 
from the color singlet field a as well. In a two-field language a non-zero a is necessarily 
generated for x 7^ by terms ax^ present in the contribution from the chiral anomaly 
for Nf = 3. 
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without having the correct asymptotic behavior 

(B) : Uoix) = HRaniSx' - W) + (1 - i?an)(x' ' '^xlx^) + Xo) (5) 

Finally, we compare our results also to a quartic potential 

(C): U,{x) = Kx' - Xlf (6) 

The quartic potential (C) is characteristic for two-flavor QCD {Nf = 2) (or 
for three flavors with small effects of the chiral anomaly). We emphasize that 
the investigation of the potentials (B) and (C) serves mainly the purpose of 
an estimate of uncertainties or errors from the unknown details of the shape 
of Uq whereas we consider (A) as a more realistic potential. 

As we have already mentioned, the dominant parameter for the character- 
istics of the chiral phase transition is Uo{0). This measures in our convention 
the potential difference between the chiral symmetric state at x = and the 
vacuum at x = Xo- It gets a contribution from the chiral anomaly f/an(0) 
and we have in our parametrization 

t/0(0) = R-JlUaniO) (7) 

Here we denote by Uanix) the contribution of the chiral anomaly, i.e. the 
part ~ Ran in eq. (^) or (|^). The size of the chiral anomaly in the vacuum, 
UaniXo) - Uan{0), is givcu for our simple potential (|) by -Uan{0). It is 
directly related |jll], |16[ to the mass of the ?7'-meson M^/ = 960 MeV and its 



decay constant fe ~ 150 MeV 



f»(0) = 45^ = ^-10-GeV* (8) 



The parameter Ran is essentially fixed for instanton-dominated multiquark 
interactions (-Ran ~ !)• This sets the scale for Uo in terms of the mass and 
decay constant of the r/'-meson. This is the basis for an interesting connection 
between the critical temperature Tc on one side and M^/ and fe on the other 
side which will become apparent later. 

The vacuum octet condensate Xo is related to the effective gauge coupling 
g and the vector meson mass flp by flp = g{fip)xo- The coupling g{fip) can 
be estimated from the p-decay width into two pions or charged leptons 
and we choose ^(770 MeV) =6. For ftp = 770 MeV this yields xo ~ 130 MeV. 
(In view of the uncertainties of this estimate one could also treat g{fip) as an 
unknown parameter within a certain range around the estimate above.) The 
coupling A is now determined by 

^ R-^f^M'^, ^ g\-^p) fjM^, 
2N]xt 2N]Ran -4 
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Another characteristic quantity for the vacuum potential is the mass of the 
a resonance given by^ 

Ml = \d'U,/dx\,. (10) 

Phenomenology indicates M^y ~ 500 MeV for realistic values of the quark 
masses. 



3 Effective masses 

We will evaluate the temperature-dependent part of the effective potential 
Af/(x, T) in a generalized mean field-type approximation. In a strong inter- 
action environment reasonable results can only be obtained if the dependence 
of the effective particle masses on the mean field and the momentum scale 
is dealt with some care. In particular, the strong running of the couplings 
should not be neglected. Furthermore, we want to be able, in principle, to 
distinguish between a first or second order phase transition. In this respect 
it is crucial that the Goldstone bosons are treated correctly, since otherwise 
a second-order transition may be obscured by an insufficient approximation. 
This is not trivial in a mean field treatment of a linear model and we devote 
some attention to this issue in sects. 9-11. 

More precisely, we include the effects of the fluctuations of the gluons, 
quarks, and pseudoscalar Goldstone bosons. Since in the low temperature 
phase gluons are associated to vector mesons and quarks to baryons [|ll 



these fluctuations include all light particles. Our approximation is formally 
equivalent to a gas of non-interacting massive particles. However, the effects 
of the interactions appear indirectly through the dependence of the particle 
masses on the mean field and the momentum scale. In particular, the gluons 
and quarks have x-dependent masses 

= /(-"k' , Ml{^,) = hl{^^)x' (11) 

where we take into account the running of the gauge and Yukawa coupling. 
The running of the couplings and the choice of the renormalization scale 
H are specified in Appendix A. From the average baryon mass in vacuum 
Mq = h^{fip)xo = 1150 MeV we extract^ h-^^i^jlp) = 9. Beyond the running 
couplings further effects of the interactions of the gluons are omitted. We 
neglect here, in particular, the thermal mass corrections which contribute 
~ gT to a chirally invariant quark mass and the gluom mass. (For simplicity, 
we also have neglected the mass splitting between the baryon octet and the 
singlet which are described by the nine quarks.) 

■^This holds if the contribution of the color singlet to the kinetic term of x is small. 
''Note that actually only Mq{x = 0) = and Mq{x = Xo) = Mq are known. The 



X-dependence of Mq{x) could be more complicated than the ansatz used in (|ll| 



8 



For the Goldstone bosons the squared mass is given for all T by the 
derivative of the potential 

..2 rri dU 36 dU 

Mh - ^aU' ^ ca— . -— (12) 

Also the effective meson decay constant depends on x 

f = c,x' = (13) 

As chiral symmetry restoration is approached for x the effective decay 
constant decreases. The same holds for the chiral condensate < ipip >~ x- 
We note that the general relation 



implies 



CGCF = — (15) 



independently of the individual constant^ cg and cp- 



4 Thermal fluctuations 

The various contributions of the thermal fluctuations are now evaluated as 

AU = AgU + AgU + AsU (16) 

with 

AgU = 2UB{gY) , \U = -UNfMhlx') , 

AM = (ATj - 1)[Jg(M2) + AJg] (17) 

Here we have defined the integrals 

J{M\T) = -{j \n{ql + M^)-j In(g2 + M2j (18) 

with /f ) = TE„/ (0, = / (0r, = (2™T)2 + g^, = ql + ^. 
For the bosonic fluctuations Jb and Jg the sum is over integer Matsubara 
frequencies n. In contrast, for Jp the Matsubara frequencies are half integer 
and therefore g|> = {2n + l)vrT, n E Z. The integrals Jb and Jp reflect the 

^The values used in (^2|) and ( p^ ) neglect the contribution of the singlet gq-condensate. 
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difference between tlie one loop expressions between nonzero and zero tem- 
perature for a given "background field" x- They are well known in thermal 
field theory 



Jb{M\T) = T j^'^ ln(l - exp(-^g2 + m^/T)) 
Jf{M\T) = T ^ ln(l + + M^T)) (19) 

and can be expressed in terms of dimensionless integrals with = M'^/T'^ 

/-oo J 

JBiF) = dqf ln(l T expi-^q^ + m^)) (20) 

We observe the Boltzmann suppression for large m, i.e. 

lim J^(^)=(+)-^ - e- (21) 



The fluctuation integral Jg for the Goldstone bosons has to be handled 
with care since Mq may be negative for a certain range of x- The Gold- 
stone boson mass Mq vanishes at the temperature-dependent minimum of 
the potential XoiT). (In our notation Xo(0) = Xo corresponds to the param- 
eter appearing in the zero-temperature effective potential (^.) The physics 
of the "nonconvex region" x < Xo{T), where our approximation leads to 
dU /d{x^) < 0, can be properly discussed |T8| only in the context of a coarse- 



grained effective action as the effective average action [T^, ^ . We give here 
some "ad hoc" prescriptions which catch the relevant physics for most (not 
all) situations. For negative U' the fluctuations with momenta < —U' 
should be replaced by spin waves, kinks or similar "tunneling" configura- 
tions which fluctuate between different minima of the potential |jl8[. Here 



we simply omit the fluctuations with q^ < —U' . As a consequence, Jg equals 
Jb for f/' > 0, whereas the lower integration boundary for q in eq. (pOD is 
replaced by \m\ if m? is negative. With 

2 j Jb{M^,T) for M2>0 

^""^^ '^^^XS^Io^ dyy^y' + \m^\Hl - e"^) for < ^^^^ 

we observe that Jg is continuous in M^. For < the contribution of 
the Goldstone fluctuations is enhanced as compared to a gas of massless free 



particles. Finally, the correction AJg in eq. ( [TTD is related to the appropriate 
treatment of the running couplings and specifled in Appendix A. 

Our treatment of the composite scalar fluctuations leads to some short- 
comings. The perhaps most apparent one concerns the behavior of AU for 
very high T. In this region the gluons have only two helicities and the factor 
24 multiplying Jb should be replaced by 16. In fact, a more correct formula- 
tion would treat the 2{N^ — 1) transversal gluon degrees of freedom separately 
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from the longitudinal degrees of freedom. The latter ones are associated to 
the scalar sector by the Higgs mechanism. In contrast to the fundamen- 
tal scalar in the electroweak sector of the standard model the octet x is a 
composite and the scalars should not be counted as independent degrees of 
freedom at high T. In a more complete treatment this is realized by large 
T-dependent mass terms for all scalars which suppress their contribution ef- 
fectively for large T. The neglection of this effect in our rough treatment 
leads to an overestimate of the bosonic contribution^ at large T by a factor 
3/2. 

On the other hand, in the low temperature region we have neglected 
bosonic contributions from the cr-resonance or the scalar octet (a^ etc.) 
which have a similar mass as the vector mesons. This results in an un- 
derestimate of the bosonic contribution which is particularly serious near a 
second-order or weak first-order phase transition. At a second-order phase 
transition all scalars contained in an appropriate linear representation be- 
come massless. In the chiral limit the minimal set are 144 real scalars for 
Nf = 3 and 32 real scalars for Nf = 2. In contrast, our treatment accounts 
only for {N^ — 1) + {Nj — 1) massless scalars at the phase transition. 

Another shortcoming is the neglection of temperature effects in the instan- 
ton-induced interactions. They become important for vrT ?a fip and can be 
neglected for low T. It is conceivable that these effects influence substantially 
the details of the phase transition. Despite of all that we remind that the 
gross features of the equation of state and the value of the critical temperature 
are relatively robust quantities. We will see below that an error of 20 % in 
the number of effective degrees of freedom in the quark-gluon phase at 
influences the value of Tc only by 5 %. A similar statement holds for a 
20 % error in the estimate of f/o(0) — f/o(xo)- For a first rough picture of the 
implications of the color octet condensate on the chiral phase transition our 
approximation of the thermal fluctuations (|T6]) seems therefore appropriate. 

5 High temperature phase transition 
for vanishing quark masses 

For a given form of the vacuum potential Uq the temperature dependence of 
the effective potential and the temperature-dependent masses of the pseudo- 
particles can now be computed. We concentrate here on a potential that is 
dominated by instanton effects reflecting the chiral anomaly and take Ran = 
0.9 in eq. (Q). In table 1 we present for the chiral limit (vanishing current- 
quark masses) the values of the vector- meson mass flp (input), the decay 
constant / and the cr-mass Mo- for the vacuum. The results for the instanton- 

^The suppression of the Goldstone boson contribution Jq at large T is included properly 
in our approximation. 
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induced potential are denoted by (A). For our parameters realistic values of 
/ip, / and Mg- are obtained for a realistic average current quark mass (see 
sect. 11, table 2). They differ from the chiral limit shown in table 1. In order 
to investigate the influence of the uncertainties in the vacuum potential on 
the thermal properties, we also compare with corresponding results for the 
potentials (B) (C) in eqs. (|^), (|^). (For the potentials (B), (C) we have not 
optimized the input parameters Jip, g{fip), hxiftp) with respect to the observed 
particle masses for nonzero current quark masses.) 

As the temperature increases, we find a first-order phase transition. The 
critical temperature Tc is indicated in table 1, as well as the Higgs con- 
tribution to the mass of the vector mesons and baryons slightly below the 
critical temperature (/i^f '^^^ (Tc) and Mj^j^^^^Tc)). We observe that the de- 
crease of the effective masses Jip(T) and Mq(T) for increasing temperature 
is largely determined by the temperature dependence of the effective gauge 
and Yukawa couplings and only to a smaller extent by the decrease of the 
expectation value Xo{T). We observe that the true screening masses get 
temperature corrections which increase ~ T and become important near Tc. 
Also the maximum of the vector meson spectral function, which is relevant 
for dilepton production, differs from the screening mass. As T increases be- 
yond Tc, the absolute minimum of U (x, T) jumps to x = and the gluons 
and quarks (or vector-mesons and baryons) become massless in our approx- 
imation. The equation of state approaches rapidly the one for a free gas of 
gluons and quarks (see sect. 6 for details). 

Since the potentials (B)(C) are qualitatively quite different from (A), we 
believe that the range Tc ~ 130-160 MeV roughly covers the uncertainty from 
the choice of the vacuum potential Uq. We note that the critical temperature 
for the instanton-induced potential (A) comes close to the value Tc = 154 ±8 
MeV suggested by a recent lattice simulation 0. If we include only the 
transversal gluon degrees of freedom for x = (see sect. 4) the value of Tc 
becomes somewhat larger by about 4 %. 



Table 1: First-order phase transition in the chiral limit {mps = 0). 



vacuum 


critical temperature 


Uo 




/ 




T 


^(^■^■^)(Tc) Mf'^'^)(Tc) 


riTc) 


A 


700 


68 


1570 


154 


290 580 


0.37 


B 


770 


113 


580 


134 


440 700 


0.66 


C 


770 


113 


480 


131 


450 720 


0.71 



6 High temperature phase 

The results of a numerical evaluation of U (x, T) can be understood qualita- 
tively by simple analytical considerations. For large temperature T > M the 



12 



integrals J are dominated by their values at M = 

Jb(0, T) = Jg(0, T) = , J^(0, T) = —T' (23) 

If the minimum of U for large T occurs at x = 0, the gauge bosons and 
quarks are effectively massless. The pressure is then given by 

We recognize the contribution of a free gas of 2{N'^ — 1) = 16 bosonic and 
ANcNf = 12Nf fermionic massless degrees of freedom. The scalar part AgU 
is subleading and will be neglected (see sect. 4). The negative contribution 
—Uq{0) = — Axo accounts in our approximation for interaction effects related 
to the octet condensation. It becomes irrelevant only for large T whereas for 
temperatures near Tc its effects cannot be neglected. For large T the pressure 
approaches rapidly the Stefan-Boltzmann limit 

7 ,^ \ OrjlA 



The energy density 

. = f/-r^ = c/„(o) + A(i + Hl^),^r' (26) 

deviates from the one for a relativistic gluon-quark gas due to Uq. One 
obtains the equation of state 

. ,rrfn^ e-3p 45 UojO) 



e + p (8 + fiVy)7r2 



This yields for Nf = 3 



r ^0192^ -r(T)^ ^-1^ f28) 
^"e + p T4 -^^^^=^4' e"3 + r ^^^^ 

We have indicated the value of t(Tc) obtained from the numerical evaluation 



of eq. (16) in the chiral limit in table 1 (and for realistic quark masses 
in table 2 in sect. 11). The small value of r for the instanton-induced 
potential (A) implies that the equation of state for a relativistic plasma 
is particularly rapidly approached as the temperature increases beyond Tc. 
The fast turnover to a relativistic plasma is consistent with findings in lattice 
gauge theories 

In a more accurate treatment the free energy of the quark gluon gas and 
the Stefan-Boltzmann law receive additional perturbative corrections also 
for large T. The dominant corrections can be incorporated in our picture 



by adding in eq. (11) the temperature dependent mass terms We 
emphasize, however, that the contribution — f/o(0) in eq. (24) remains as an 
important non-perturbative effect. 
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7 Chiral phase transition 



For a first-order transition the critical temperature Tc can be related to Uo{0). 
Let us denote the value of the potential in the hadronic phase at Tc by 
Ussb = U{xo{Tc))- Equating this with eq. ( P^D yields 



= TT^ 1 + ^ mO) - Ussb) (29) 



^7r2 ^ ^ 32 

where we have neglected A^f/ since the Goldstone boson fluctuations are 
not relevant at x = for a sufficiently strong first-order transition. As long 
as Tc remains much smaller than the mass of the gauge bosons ft^^^^^^ and 
fermions M^^^^^ in the phase with broken chiral symmetry, their contribution 
to Ussb remains strongly Boltzmann-suppressed (cf. table 1). For a strong 
first-order phase transition the quantity Ussb is therefore dominated by the 
fluctuations of the Nj — 1 Goldstone bosons. If Xo(^) is close to Xo = Xo(0) 
we can also neglect the difference Uo(xo{T)) — f/o(xo) and Ussb is given by 
the contribution of a free gas of Goldstone bosons 

Ussb = (30) 

As a consequence, the critical temperature is given by Uo{0) quite indepen- 
dently of the details of the shape of Uo{x) or the thermal fluctuations 



Ti = ^(l + ^-^)-^f/o(0) = 0.231 f/o(0) (31) 



(For the quantitative estimate we use Nf = 3.) For the approximation (^) 
this yields an equation of state 

— ^ ^ r(Tc)^ ^ 0.83^ (32) 

and we see again that a free relativistic gas of gluons and quarks is reached 
rapidly as T increases beyond T^. The pressure is continuous at Tc and 
directly given by p(Tc) = —Ussb- Comparing eq. (|30|) with the Stefan- 
Boltzmann value (H) yields p{T^)/psB = § (^) for Nf = 3(2). 

Eq. ( ^11) relates the critical temperature to the properties of the ?7'-meson 
and the fractional contribution Ran of the chiral anomaly to Uo{0) (cf. eq. 
(i) 

Tc = R-y^NJ^^^ ■ 221 MeV (33) 

For Ran = (0.3, 0.5, 0.8, 1) and Nf = 3 this yields Tc = (173, 152, 135, 128) 
MeV. Fluctuations in the SSB-phase beyond the pions (30) further lower the 
value of Ussb- This replaces in eq. (33) Ran by R, with R < Ran, leading to 
an increase of T^. For three-flavor QCD our numerical results are consistent 
with this qualitative picture of a strong first-order phase transition. 
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The quantitative estimates (pl|)- (|3^) are only valid if the minimum at 
Xo(Tc) is still near xo- This depends on the influence of the Goldstone bosons 
and, in particular, on the ratio between and the meson decay constant / 
(see sect. 11). For a weak first-order transition or a second-order transition 
(for Nf = 2) the potential difference f/o(Xo(^c)) has to be included in Ussb 
(PPI). This contribution lowers the critical temperature. On the other hand, 
the gauge boson and quark fiuctuations may not be negligible anymore for low 
enough xo{Tc)- This effect enhances the critical temperature. The difference 
between the rough estimates (|33|) and (^) and the numerical conputation 
(table 1) is related to these effects[|. It is particularly pronounced for the 
potential (A) for which the gauge boson mass in the SSB-phase is already 
considerably smaller at than at T = 0. 

In the two-flavor case the chiral phase transition in QCD ressembles in 
many respects the electroweak phase transition. The vacuum potential is 
analytic in x^ at = such that a flrst-order transition can only be induced 
by the fluctuations of the gauge bosons. In this analogy the octet replaces the 
Higgs-scalar and plays the role of the mass of the Higgs scalar Mh- Also 
the gluons replace the weak gauge bosons and the quarks play the role of the 
top quark in the electroweak theory. The electroweak phase diagram shows 
a flrst-order phase transition for small Mh and one may expect the same 
for QCD for small M„. In the electroweak theory this transition ends for 
some critical Mh,c in a second-order transition, whereas for Mh > Mh,c the 
transition is replaced by a continuous crossover 0, ||2^,p3[. In the case of 



QCD with vanishing quark masses one has, however, an additional important 
ingredient, namely the existence of an order parameter and Goldstone bosons. 
In fact, in the absence of current quark masses an order parameter Xo(^) 7^ 
breaks spontaneously a global symmetry. In consequence, there must be a 
phase transition as Xo(^) reaches zero when the critical temperature is 
reached from below. This holds for arbitrarily large Mg.. In two-flavor QCD 
with vanishing quark masses a second-order phase transition replaces the 
crossover of the electroweak theory. For nonzero quark mass the analogy 
is even closer, with a change from a flrst-order transition to crossover in 
dependence on M^. In contrast to the electroweak theory, however, M^- is not 
a free parameter but can, in principle, be computed in QCD. An important 
ingredient for the determination of the characteristics of the phase transition 
is the size of the gauge coupling which is large in QCD. The realization of 
a second-order transition/crossover seems therefore quite plausible in two- 
flavor QCD. 



''The numerical computation also uses 24 gauge boson degrees of freedom (eq. (|l7|)) 
instead of 16 for the analytical discussion in sects. 7 and 8. This lowers Tc by 4 %. 
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8 Spontaneous chiral symmetry breaking for 

T <T, 

For an understanding of the temperature effects in the hadronic phase for 
T < Tc we need the value of Xo(T). The octet expectation value in thermal 
equilibrium obeys the field equation 

dU dlJ 

^x-'-w^r'-"'-'^ (34) 

where = in absence of quark masses. For a determination of the temper- 
ature dependence of the order parameter Xo(^) it is useful to investigate the 
field equation analytically. For the gluon contribution one finds from (|T7|) 

A,f/' = Q^^^oU = l2g\^,,{T))lf\^,l{nT) (35) 

Here 

reflects the running of the gauge couphng in dependence on the temperature- 
dependent gluon mass 

f^%T) = g'if^Dx' (37) 

The integral 

/f)(M^T) = ^Q^MM',T) 

d?q 1 f d^q 1 



(27r)3 (f -f (27rnT)2 + M2 J {2^^ + 



27r2 Jq " ey+'^ - 1 



contains for large T terms quadratic, linear and logarithmic in T 

12 47r + 167r2 ^^^^ 

More precisely, the expression ( |39D applies for ^ vr^T^ and we will 
drop the logarithm in the following. We observe an infrared divergence in 
dll^^/dM^ for M 0. 

The fermionic contribution reads 

^ A,U = -6Njh'll''\hlx',T) (40) 
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with 



= ^(^x(/^t)x^) (41) 



The corresponding integral 



T2 /•°o^^yy^ + 2my 1^2 



, , dy^^^^ = —T' + ... (42) 

involves half-integer Matsubara frequencies and does not show an infrared 
divergence. 

The Goldstone-boson integral if^ = 2dJG/dM'^ equals If for > 0. 
For negative it becomes 

/r = -^y^ c^l/Z/(l/'+|m^l)^'/'ln(l-e-^) = - + Ajf) (43) 
For small values of \M\/T one finds 

.,.-^^(.n^-l).... (44, 

and we see that jj*^^ is continuous at = 0. On the other hand, for large 
\M\/T the resuh is 

/P^ = + ... , c = /°° dy-^ ^ 2.4041 (45) 

^ 47r2|M| ' io e^-l ^ ^ 

For sufficiently large T a positive temperature-dependent mass-like term 
dominates the potential derivative at the origin 

^ =(e + ^hAT- (46) 

This overwhelms any classical contribution and there is therefore no nontriv- 
ial minimum of U{x,T) for x 7^ 0. Chiral symmetry is restored {xo{T) = 0) 
and the gluons do no longer acquire a mass from the octet condensate. On 
the other hand, for T < Tc the absolute minimum of the potential occurs for 
Xo(T) 7^ 0, corresponding to a nontrivial solution of eq. (|3^). 
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9 Goldstone-boson fluctuations 



For a determination of Xo{T) in the low temperature phase with chiral sym- 
metry breaking the Goldstone boson fluctuations are important. They dom- 
inate the low-temperature behavior and play an important role at a possi- 
ble second-order phase transition. The correct treatment of the Goldstone 
bosons in a mean field-type approach is rather subtle due to complications in 
the infrared physics for massless bosons. (This also applies to gauge bosons.) 
We show in the next section that our renormalization-group motivated pre- 
scriptions lead to the correct behavior at low T and we demonstrate in ap- 
pendix B that they also can account^ for the possibility of a second-order 
phase transition. The latter is particularly relevant for the case of two light 
quark flavors or for realistic QCD where the physical strange quark mass 
seems to correspond to a very weak first order transition or a crossover, in 
the vicinity of a second order transition for an nearby critical strange quark 
mass. 

The Goldstone boson mass involves the derivative of the temperature- 
dependent effective potential Mq = cqU', which obeys for all T 

U' = B{x) + ^^^ccU"/\lf\caU\T) + ... (47) 
B{x) = U', + -^ccT^U'^ + Wl'^\g\\ T) - 6iV,/^2/^(/^Jx^ T) 

Here the dots stand for corrections arising from A Jg (|A.3|) which vanish for 
f/' = and are negligible in the vicinity of an extremum. /^From Aj}*^^ (0, T) = 
one infers a simple equation for the extrema of U away from the origin. 
The condition [/' = is realized for 5 = or 

where we recall the definitions 



(48) 



U' - U" - ^ - (49) 

'~2xdx ' ""'-AxAdx' xdx ^ ^ 



For low temperature the r.h.s. of eq. (^81) is Boltzmann-suppressed and can 
be neglected. 

The differential equation (|47| ) for U'{x) can be turned into an algebraic 
"Schwinger- Dyson equation" by replacing U by t/g . This is reasonable for 
Uq > 0, whereas for negative Uq we omit the contribution ~ Uq in B and 

® A naive mean field treatment typically produces a first-order transition when massless 
bosons are present. This is due to the nonanalytic "cubic term" in the expansion of Jb 
(eq. ( p^ ) in powers of M. 
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in eq. ( ^7|) such that U' = B. These prescription^ determine U' for 
given values of x ^"^^ T. Not too far from the minimum one may further 
approximate (for Uq > 0) 

U' = B- ^^^cfu;;T ^^Wem - l^-U'&{-U'){\n{2\UVT) - 1)) 

(50) 

A reasonable approximate formula for U' used for our numerical computation 
is given by 

U' = ^sign{B){^E^ + A\B\- Ef, 

E = ^^cf Tf/^'e(f7^') (51) 

This implies that for small B in the vicinity of the minimum where B ~ 
X - Xo, one hasQ t/' ~ (x - Xof- 



10 Low temperature pion gas and chiral per- 
turbation theory 

The thermodjTiamics for low temperature should be described by a gas of 
interacting pions. Any analytical computation which pretends to describe the 
transition from a pion gas at low T to a quark gluon plasma at high T should 
reproduce the low temperature limit correctly. Indeed, in our formulation 
only the massless or very light particles contribute for low T. In particular, 
around the minimum of U near xo both the baryons and the vector mesons 
(or, equivalently, the quarks and the gluons) are heavy and exponentially 
Boltzmann-suppressed. Only the pions are light and one expects that they 
completely dominate the temperature effects at low T for x near xo- 

The effects of a thermal interacting pion gas are described by chiral per- 
turbation theory |25| which predicts for the chiral condensatej^ 



XojT) ^ <^^> (T) ^ _ Nf-1 _ N]-l (T^y 

Xo(0) <^^>{0) Nj 12/2 2N] [up) ^ ^ 

Within the nonlinear setting of chiral perturbation theory one can understand 
the temperature effects of the interacting pions on the chiral condensate in a 



^In regions of the potential where the gauge boson and fermion fluctuations are im- 
portant a better approximation would replace Uq Uq + AgU" + A^C/", similar to eq. 



J)- 

^°This holds for E > 0. Note that for U' ^ x ^ Xo the nonanalyticity in A/-[*^-' would 
destabilize the minimum. 

^"'^Here we have neglected the T-dependence of the wave function rcnormalization which 
relates x find < ipip >. 
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staightforward way. In particular, it is obvious that the lowest two orders can 
only depend on the ratio T// In order to describe the phase transition 
and the high temperature phase we are bound, however, to use a linear 
description for the scalar fields. The correct reproduction of lowest order 
chiral perturbation theory is not trivial in a linear setting and can be used 
as a test for the correctness of the particular formulation of the meanfield 
approximation. 

In our picture we can parametrize the effective potential in the vicinity 
of the minimum xo by {U = Uq + AU) 

Uo = ^M'Ax-xoY 

AU = (iVj-l)(-^T^ + gTWo/9(x^)) (53) 
Up to the x-independent contribution ~ this yields 

U = ^M^Aix - Xof + '^ ^(X - xo)} (54) 

and one infers that the T-dependent minimum is independent of 

XojT) CgCf. 2 .^ Xo(0) 

^ = 1 - — ev/ - 1)1275^ 

With cgCf = 4:/Nf (cf. eq. (|T3])) we can verify that the lowest order of an 
expansion in (xo(r) — Xo(0))/xo(0) corresponds indeed to chiral perturbation 
theory. This shows the consistency of our treatment of the scalar fluctuations 
as discussed in Appendix A. A lack of care in the treatment of the scalar 
fluctuations easily leads to inconsistency with chiral perturbation theory even 
in lowest order. In fact, it is crucial that the term ~ in AU (eq. (^) 
involves the derivative of Uq without additional temperature corrections to 
the potential. 



11 Nonzero quark masses 

We flnally discuss the effect of the quark masses in the limit where they are 
all equal. A quark mass term adds to the potential a lineaiQ piece 

Uj = Uo + Um , Um = -Jxix - Xo) (56) 

^^We omit here nonlinear quark mass corrections from the chiral anomaly. We also 
observe that the linear term actually appears for the gg-color singlet a and is transmuted 
to X only by integrating out a. This results again in a nonlinearity of Um- We consider Um 
as an approximation for T <Tc. There is actually only a singlet and no octet condensate 
in the high temperature phase for nonvanishing quark masses. 
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Therefore the location xo of the minimum oiU + Um obeys for arbitrary T 

^(Xo)=Jx (57) 

Here is proportional to the current quark mass and can be related to the 
mass of the octet of light pseudo-Goldstone bosons (tt, rj) by 

rr^PsiT) = CG^(Xo(T)) = (58) 

In particular, xo — Xo(0)is determined by the vacuum mass nips — mps{0) 

%{Xo) = mls/cG (59) 

Vanishing pressure in the vacuum (at T = 0) requires C/o(Xo) + t4i(Xo) — 
or 

Xo = Xo - Uo{xo)/jx (60) 
For a small quark mass (small j^) both xo and xo are given approximately 
by Xo- 

Near the chiral limit we can understand various quark mass effects on the 
chiral phase transition analytically: 

(1) The vacuum expectation vahic xo increases and Uo{0)+UmiO) acquires 
an additional positive contribution j^Xo- Both effects enhance Tc. 

(2) In the high-tcmpcraturc phase the minimum occurs at Xs{T) 7^ 0. 
For a first-order transition this reduces the difference between Uj in the 
low and high temperature phases and therefore diminishes somewhat the 
increase of Tf.. Also the quarks get a nonzero mass in the high-temperature 
phase, h^Xs{T)- This decreases the pressure of the fluctuations at given T 
and therefore enhances the value of needed in order to compensate the 
difference in Uj . For small this effect is only quadratic in j-^ and the linear 
effect in Uj{xs{T)) - Uj{xo{T)) dominates. 

(3) As Tc and Xs(^c) increase with increasing j^, the difference between 
the two local minima characteristic for the flrst-order transition becomes less 
and less pronounced. Finally, the transition line ends for a critical j^^c at a 
second-order transition, with crossover for > j^^c- The situation is similar 
to the liquid-gas transition. 

(4) In case of a second-order transition for — the chiral transition 
turns to a crossover for all nonzero j^. 

In summary, nonvanishing quark masses enhance the critical temperature 
and make transitions less pronounced. The case of realistic QCD with rris 7^ 
rriu^d is somewhat more complicated since the expectation values in the high 
temperature phase differ in the strange and nonstrange directions. In case 
of a strong flrst-order transition this modiflcation is, however, only of little 
quantitative relevance. On the other hand, one may envisage a situation 
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where for 7^ 0, m„ = = the critical behavior ressembles two-flavor 
QCD. 

Our numerical evaluation of U{x, T) is easily exploited for nonvanishing 
quark masses. They only enter into the determination of the location of 
the minima at Xs{T) and Xo(T) via eq. (|58D . We indicate the vacuum 
properties for anomaly-dominated potentials {Ran = 0.9) in table 2 for the 
same values of the input parameters as used in table 1. Comparison between 
the two tables provides an estimate of the effect of nonzero degenerate quark 
masses. We employ here a "realistic" average mass = (390 MeV)^ ~ 
|(2M|- -|- M^) which corresponds to the neglection of S'f/(3)-violating mass 
splittings in the pseudoscalar octet. We observe a sizeable increase of / = 
\j7/9xo and a corresponding increase in the baryon mass0 Mg = h^{fip)xQ 
and vector meson mass flp = g{fip)xo in the vacuum. 

For all three potentials (A)(B)(C) we find a first-order transition for 
mp5=390 MeV. The critical temperature is increased substantially as 
compared to the chiral limit, as can be inferred from table 2. 



Table 2: First-order phase transition for realistic average quark masses 

(mps = 390 MeV). 



vacuum 


critical temperature 


Uo 




/ 




Tc 


fi^p^'^KTc) Mf'^'«)(T,) 




A 


770 


116 


460 


170 


290 600 


0.53 


B 


800 


132 


770 


180 


470 800 


0.77 


C 


810 


142 


660 


183 


490 840 


0.77 



The mass effect on the critical temperature is moderate for our instanton- 
induced potential (A) for which the increase of about 10 % is roughly consis- 
tent with a recent lattice simulation |^. In this respect the instanton-induced 
potential computed in fl^ does apparently much better than the two polyno- 
mial potentials (B)(C) that we use for comparison. We point out, however, 
that the neglection of the quark mass effects for the instanton-induced inter- 
action leads to a substantial uncertainty in the quark mass dependence at 
the present stage. 



12 Z3- Vortices 

In an attempt to find distinctive features of our scenario which could be tested 
by lattice simultations, we may look at topologically stable excitations. The 
octet condensate in the low temperature phase indeed implies the stability 

^■^Our choice of the Yukawa coupling leads to too large values of the baryon mass in 
vacuum for the potentials (B)(C). This is only of little quantitative relevance for the results 
presented in this note and can be corrected easily by the choice of a smaller Yukawa 
coupling. 
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of macroscopic Zs-vortices. They are absent in the high temperature phase. 
In the early universe, such strings would have been produced as topological 
defects during the QCD-phase transition. The string tension is typically 
of the order a ~ such that Ga ~ {mp/MpY is tiny. Observational 
consequences of the production of typically one string per horizon at the 
time of the phase transition are not obvious - in particular such strings are 
not candidates for seeds of a later galaxy formation. 

As topological defects, the Za-vortices correspond to the nontrivial ho- 
motopy group Tiii^SUi^) / Z^) = Z3. The color octet (and singlet) scalars 
transform trivially under the center of the color group, similar to the glu- 
ons. Any classical bosonic field configuration is therefore invariant under 
Za-transformations. Once the color group gets "spontaneously broken" by 
the octet condensate, the Za-vortices become topologically stable. As an 
example, consider a static vortex in the z-direction with octet scalar field 

Xijab = x{r) iXz)ab (v\^)X,v{<^)),, (61) 

Here v is given by a y^-dependent gauge transformation 

v{<^) = exp[^<^Xs) (62) 

such that a rotation around 2tt corresponds to an element of 

,27rz\ , /27ri, . 

expi-^Xs) = exp{—) (63) 

The scalar field is therefore well defined and free of singularities provided 
x{r) behaves properly at the origin, e. g. %('" = 0) = 0. For r — 00 we 
assume that x{^) approaches the expectation value xo which minimizes the 
effective potential, xi"^ ~^ 00) Xo- The potential energy of the string 
is therefore concentrated in the core of the string, typically of radial size 
rrip. The nontrivial homotopy group tells us that no smooth deformation can 
change the nontrivial behaviour for r 00. In particular, no (singularity 
free) gauge transformation can "unwind" the string. 

Without the gauge fields, the gradient energy per length of the configura- 
tion (0), namely Egrad ~ I drr-'^d^XijabdipXijab ~ J drr-'^x'^{r), would still 
diverge logarithmically for r — 00. We may supplement a gauge field in the 
(y9-direction which becomes a pure gauge for r — > oo(^a{r — > 00) = l,a(r = 

0) = 0) 

i 1 

= — d^v'^ v*a{r) = Xsa{r) (64) 
9 V3g 

It is then easy to see that the covariant derivative vanishes for r — ^ 00 

D^Xijab = d^Xijab - 'i9{Ak)vXkjab + WXikab^Akj) ^ ^ 0. (65) 
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The gauge field strength vanishes as well in this limit. We are therefore sure 
that the vortex has a finite string tension, the precise value depending on 
the shape of the functions x(^) and a(r) which have to be determined by 
solving the field equations with the appropriate boundary values for r = 
and r —>■ oo. 

Actually, the topological situation which leads to stable vortices gets 
somewhat more involved by the fact that both color and flavor groups are 
broken simultaneously by the octet condensate. A color rotation in the 
element of SU{3)c can be "unwound" by a yj-dependent color-fiavor-locked 
rotation. This shifts the nontrivial topology from the color sector to the 
flavor sector, more precisely to the nontrivial center of the diagonal flavor 
transformations from SU{3)l x SU{3)fi. The topological stability due to the 
nontrivial homotopy group tti is not affected. Since QCD is invariant only 
under global flavor transformations the vortex is, however, not invariant un- 
der coordinate-dependent color-flavor-locked transformations. The vanishing 
of the gauge covariant kinetic term for r ^ oo occurs only if the 93-dependent 
rotation of the vacuum state is associated to the color sector. 

In this context it is useful to recall the global structure of the symmetry 
group of QCD with three massless quarks. For infinitesimal transformations 
the symmetry group is 

G = SU{3), X SU{3)l X SU{3)r x U{1)b (66) 

Global rotations in the center of SU (3)c and the center of the diagonal vector- 
like flavor group SU{3)v correspond to identical phase rotations of the quark 
fields. These transformations belong also to the abelian group U{1)b which is 
associated to conserved baryon number. A direct consequence of this global 
group structure is the triality rule: Color representations in the class (3, 6, 
etc.) have baryon number B = 1/3 mod 1 and electric charge Q = 2/3 mod 1 
whereas the class (1, 8, etc) has B = 1 mod 1 and Q = 1 mod 1. This ex- 
plains why the physical fermions are baryons with B = 1 and integer Q, since 
all physical states are color singlets. (Some care is needed for the correct in- 
terpretation of the Higgs picture |TT|.) 

In a first look it may seem that the macroscopic Z3- vortices are not present 
in the confinement picture of QCD, in contradiction to our assumption of 
equivalence of the Higgs and confinement pictures. This is, however, not 
obvious. In gluodynamics (QCD without light quarks) the ZAr-vortices are 
discussed as important configurations for the understanding of confinement 



2£] . Lattice simulations could establish the existence of macroscopic vortices 
in the vacuum of three flavor QCD and their absence in the high temperature 
phase. Such a finding (or the contrary) may be interpreted as an important 
test for our picture. 
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13 Conclusion 



c 



In conclusion, we have presented here a rather simple picture of the high- 
temperature phase transition in QCD. Both confinement and chiral symmetry 
restoration are associated to the melting of a color-octet condensate at the 
critical temperature. The main phenomenological features of our picture of 
the phase transition are the following: 

(1) For three light quarks with equal mass a first-order phase transi- 
tion separates a low temperature "hadronic phase" from a high temperature 
"quark-gluon phase". In the chiral limit of vanishing current quark masses 
the quarks and gluons are massless in the high temperature phase where 
chiral symmetry is restored. Color symmetry becomes a symmetry of the 
spectrum of pseudoparticles above the critical temperature Tc. Below 
both chiral symmetry and color symmetry are spontaneously broken. Chi- 
ral symmetry breaking generates a mass for the quarks which now appear 
as baryons. Eight massless Goldstone bosons signal the global symmetry 
breaking. Local color symmetry breaking gives a mass to the gluons by the 
Higgs mechanism. The gluons are identified with the octet of vector mesons 
{p,K*,uj). The first-order transition is therefore associated with a jump in 
the vector meson mass (they become massless gluons for T > Tc) and the 
baryon mass (they turn to massless quarks for T > Tc) . Our picture of gluon- 
meson duality and quark-baryon duality allows us to describe the excitations 
relevant above and below Tc by the same fields. The first-order phase tran- 
sition extends to nonzero current quark mass, including values which lead 
in the vacuum to a "realistic" average mass for the (tt, K, ?])-pseudoscalars 



mps = 390 MeV « •^(2M|- + M^)/3. 

(2) The effective multiquark interactions at low momentum or the corre- 
sponding effective potential for gg-bilinears are presumably dominated by the 
axial anomaly arising from instanton effects. We can then relate the critical 
temperature of a strong first-order transition to the mass and decay constant 
of the meson. For mps = 390 MeV we find Tc ~ 170 MeV whereas the 
critical temperature is lower in the chiral limit. For an anomaly dominated 
vacuum potential the prediction for the chiral limit {mps = 0) is independent 
of many details and found in the range Tc ~ 130 — 160 MeV. 

(3) A dynamical picture for an instanton-induced color octet condensate 
has been developed recently [0 . This has led to a computation of the octet 



potential which involves the running gauge coupling as the only free param- 
eter. The characteristic features of the instanton-induced vacuum effective 
potential are depicted by the potential (A), eq. (H). We concentrate in 
the following on this estimate. It is encouraging that the predictions of the 
instanton-induced potential for the critical temperature, namely Tc = 170 
MeV for mps = 390 MeV and Tc = 155 MeV for mps = 0, agree well with 
recent lattice simulations [0. One may interpret this as a test for the idea of 
spontaneous color symmetry breaking by instanton effects. 
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(4) For the equation of state in the quark-gluon phase we find at Tc a 
ratio between pressure and energy density p/e ~ 0.13 for mps = 390 MeV 
and p/e ~ 0.19 for rups = 0. This ratio approaches very rapidly the equation 
of state of a relativistic gas {p/e = 1/3) as T increases beyond Tc. 

(5) In the hadronic phase the screening masses of baryons, vector mesons 
and pseudoscalars show a strong temperature dependence as T approaches 
Tc from below. Near Tc the chiral symmetry breaking contribution to the 
average baryon mass Mg ^ 600 MeV is only about one half of the vacuum 
mass and approaches a typical constituent quark mass. Also the Higgs con- 
tribution to the p-meson mass is found substantially smaller than in the 
vacuum, fip{Tc) ~ 300 MeV. For T > Tc the effective chiral symmetry break- 
ing fermion mass jumps to a small value Mg ^ 40 MeV which is compatible 
for rups = 390 MeV with common estimates for the average current quark 
mass rriq = ^{rris + rrid + The average pseudoscalar mass is found as 
mps{Tc) ~ 180 MeV. It is therefore substantially smaller than in the vacuum. 

The present paper also contains a rather detailed discussion of a mean 
field-type computation of the temperature dependence of the effective po- 
tential. This has become necessary since a too naive treatment of the scalar 
fluctuations has previously often led to incorrect results near second-order 
or weak flrst-order phase transitions. The reason is the complicated infrared 
behavior of massless boson fluctuations in thermal equilibrium. A second 
problem is the non-convexity of the perturbative effective potential which 
leads naively to negative scalar mass terms. Most previous work has simply 
left out the scalar fluctuations because of the technical difficulties associated 
with these problems. Nevertheless, the Goldstone boson fluctuations are 
important for QCD at low temperature and near the phase transition, and 
we therefore want to include them. Our prescriptions are based on lessons 
from earlier more involved renormalizat ion-group investigations and should 
include the most dominant effects. In particular, they account correctly for 
the low temperature behavior as described by chiral perturbation theory and 
are consistent with the physics of a second-order phase transition if this oc- 
curs. Our treatment can easily be taken over for other models. 

Our approach to the high-temperature phase transition in QCD still con- 
tains many uncertainties, which have been discussed at various places in this 
paper. In particular, a better understanding of the effective potential in the 
vacuum and a separate treatment of the color octet and singlet condensates 
would be most welcome and probably needed for a study of flner aspects of 
the phase transition. This concerns, in particular, the interesting question 
about the existence and order of a phase transition in "real QCD" with differ- 
ent strange quark and up/down quark masses. An appropriate tool for this 
purpose seems to be a renormalizat ion group study similar to the succesful 
treatment of the phase transition in the Nambu-Jona-Lasinio model. 

Nevertheless, we flnd the simple overall picture of the chiral and decon- 
flnement transition rather convincing and the quantitative results of a flrst 
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rough computation encouraging. We hope that this work can become a 
starting point for a quantitative analytical understanding of the QCD phase 
transition. 



Appendix A: Running couplings 



In general, the effective masses M appearing in the J-integrals ( |18D depend 
on the momentum (f. This could be expressed through a momentum de- 
pendence of the effective gauge and Yukawa couplings. We want to avoid 
here the complications of momentum-dependent masses. Instead, the dom- 
inant effects of the running couplings are taken into account by the choice 
of an appropriate renormalization scale /x in eq. (jl^). Corrections to this 
approximation are reflected in the term AJc in eq. (p!?!). 

For not too large m the integrals (|I^ are dominated by momenta ~ 
(vrT)^. We therefore may choose an appropriate temperature-dependent 
renormalization scale 



(Here //^ = /U^(/i = fj,p) refers to the vacuum with arbitrary x corresponding 
to a possible presence of "sources" or quark mass terms.) This choice of ht 
is a valid approximation for the fermion fluctuations and we therefore use 
h'^{fiT) in the argument of Jp in eq. ([T7|) . 

For bosons the issue is more involved due to the infrared behavior of 
classical statistics. The first two terms in a Taylor expansion, Jb(0) and 
9Jb/(9M|^j2=q, are dominated again by momenta (f ~ {tiTY. The remain- 
ing integral Jb = Jb — Jb{0) — ^^f^|^j2=o receives, however, an essential 
contribution from the n = Matsubara frequency which corresponds to 
classical statistics. The three-dimensional momentum integrals of classcial 
statistics are more infrared- singular than the four-dimensional integrals for 
the vacuum. In consequence, the remaining integral Jb is dominated by mo- 
menta (f ~ M^. This leads in perturbation theory to the so-called cubic 
term |^ which has played an important role in the discussion of the elec- 
troweak phase transition p8|. Since Jb is not analytic at = 0, a careful 



treatment of the long-distance physics is mandatory if one aims for quantita- 
tive precision. For first-order phase transitions the momentum dependence 
of has only moderate effects for the gauge boson fluctuations^^ and we 
choose = 5'^(/ir)x^ in Jb- 

For the Goldstone boson fluctuations more care is needed if one wants 
to be consistent with chiral perturbation theory for low T and the correct 

""^■^For a second-order phase transition one may split the momentum integral in Jb and 
use = g^{iJ,p)x^ for (j^ < (a^p/2)^- Here grifJ-p) obeys an effective three-dimensional 
evolution for /Xp < ttT . 
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behavior at a second-order phase transition. We first include the gauge bo- 
son and fermion thermal fluctuations. At this level the mass term for the 
Goldstone bosons is given by 

^G,o = ^^^(^0 + A,t/ + A,U) (A.2) 

In a second step we take the scalar fluctuations into account. Then only 
the complete mass term Mq (|T2D vanishes at the minimum of U{x,T). The 
thermally corrected (pseudo) particle mass term Mq is relevant for momenta 
(f ~ \Mq\ and we use Mq for the argument of Jg- On the other hand, the 
high momentum part of the integral Jq is dominated by momenta where the 
thermal mass corrections from the scalar fluctuations are not yet important. 
Therefore Mqq ([A.2|) is relevant in this momentum range. We account for 
this by the correction 

^Jg = ^{M'q,o - Ml)fQ{Ml/T') (A.3) 

which effectively replaces Mq by Mq q for the term ~ which is dominated 
by momenta (f ~ (ttT)^ if > \Mq\. On the other hand, the factor fc 
reflects the fact that all temperature fluctuations and therefore also A are 
exponentially Boltzmann-suppressed for ^ Mq. We choose here a simple 
form (7 = 0.5) which mimicks eq. ( pi]) for large im? 

fcim') = exp(7 - (7^ + (m^)^)^/^)!! + ^^y^^ (A.4) 

T 

These prescriptions may seem somewhat ad hoc and technical. They find a 
deeper motivation from the study of the renormalization flow in thermal field 
theories |ll9|, pO|. We will see below that they guarantee agreement with chiral 
perturbation theory for low temperature and avoid unphysical singularities. 

We finally have to specify the scale dependence of the gauge and Yukawa 
coupling. For the running gauge coupling we employ the perturbative 
function in three-loop order in the MS scheme 

' l3o = ll~2Nf/3 (A.5) 

Our normalization g{770 MeV)=6 is suggested by the phenomenology of 
p-decays into pions and For the Yukawa coupling we use the 



one-loop renormalization-group equation 
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Treating also the gauge coupling in one-loop order this has the solution 



ahl 



(2/3o - h)g\^i,) 



2/3o' 



-1/2 



(A.7) 

We use the approximation ([A.7|) with b = 16 such that we recover for small 
the standard anomalous dimension for the mass and take somewhat arbitrar- 



ily]^ a = 1. The normalization is fixed by eq. 
GeV/xo- 



0), i.e. (770 MeV)=1.15 



Appendix B: Second-order phase transition 

For an investigation of a possible second-order phase transition in the chiral 
limit we concentrate on the behavior of U near = or small values of 
X^/T'^- In the vicinity of the critical temperature of a second-order (or weak 
first-order) transition the Goldstone fluctuations play a role in this region. 
The derivative of the temperature-dependent effective potential becomes (in 
the range where U' > 0) 



+ ^S^c^f/Vt/' 1 - (B.l) 




Here we have improved our treatment of the gauge boson fluctuations by us- 
ing g{^p) instead of gifJ^r) for the low momentum fluctuations corresponding 
to classical statistics (the term ~ T). Let us denote by Tq the tempera- 
ture where U' vanishes at x = 0, i.e. U'{0,Tq) = 0. A second-order phase 
transition at = Tq occurs if f/'(x, To) is strictly positive for all x > 0. 
In contrast, one has a first-order transition if [/'(x, Tq) is negative for small 
X > 0. Then the point x = corresponds to a maximum of U{x, To) and the 
critical temperature for the first-order transition is below Tq. 

Consider first the potential (C) which is relevant for the two-flavor case. 
Omitting for a moment the Goldstone boson contributions and the running 
of the gauge coupling {Pg = 0), the term ~ g^xT in eq. ( p.l|) would lead to a 
first-order transition. This corresponds to the well-known "cubic" term in the 
effective potential for the Higgs scalar in the electroweak theory. As has been 
discussed extensively in connection with the electroweak phase transition 



2^, the infrared physics of the gauge bosons responsible for this term has 
to be handled with care. In particular, the running of the gauge coupling 
is crucial for a correct picture. For a running g the product gx = /ip goes 

^^A more accurate treatment would keep the color octet and singlet as separate fields. 
The evolution equations for the two respective Yukawa couplings differ. 
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to a constant A(T) for x ^ and the nonanalyticity in disappears. We 
note that x-independent terms in eq. ( [B.l| ) only influence the location of Tq. 
Denoting f/"(0, Tq) = one obtains for small x ^i-nd \T — Tq\ an expansion, 

-2 



with constants Ci, C2 and A^? - 1 = iV 



Nc 



3/2 



t/' = gxX' + Ci{T' - T2) + C2(T - To) - — i^^^TVt/' (B.2) 

OTT 



For T ^ To and x ^ this yields ||26 

Stt 



t/'(p,T) 



^3/2^ (x^ + C3(T-To)) 
^cj To 



(B.3) 



and describes |^ a second-order phase transition. The critical exponent u = 1 
correspond^ to the leading order of the expansion in the effective three- 
dimensional theory |2^. An improvement of the description of the critical 



behavior (with more realistic critical exponents) needs to incorporate the 
fluctuations of the neglected scalar modes (see sect. 4) and can be done with 
the help of modern renormalization group methods []T^. We conclude that 
a second-order chiral phase transition for Nf = 2 is compatible with our 
picture, provided the infrared fluctuations of the gauge bosons are treated in 
the appropriate way. 

For three massless flavors the presence of an effective cubic term ~ x^ 
in the effective potential at zero temperature (eq. (|^)) changes the situation 
profoundly. This term reflects the axial anomaly and is induced by instanton 
effects. The term ~ — x always dominates the r.h.s. of eq. ( |B.1| ) for small 
enough values of x since there is no competing term linear in x- One infers 
the existence of a first-order phase transition for three massless flavors. 

The vicinity of the critical temperature of a second-order or weak first- 
order transition is governed by universal behavior. The universality class 
is characterized by the symmetry and the representations which remain (al- 
most) massless at T^. For two-flavor QCD it is not obvious which are the 
massless scalar representations at a second-order phase transition. Further- 
more, the gauge bosons behave essentially as a free gas near the transition. 
Even though their low momentum classical fluctuations may acquire an ef- 
fective mass, they will influence the nonuniversal critical behavior. This 
influence of the gauge bosons is characteristic for a simultaneous chiral and 
deconfinement transition and leads to modifications of results obtained in 
quark- meson or Nambu-Jona-Lasinio models [^. In view of the distance 
scales probed in lattice simulations, we find it rather unlikely that critical 
behavior in the standard 0(4)-universality class will show up there. This is 

^^A first-order transition is possible if a new minimum appears before the running g has 
reached the behavior gx = const.. 

"'^''This differs from high-temperature perturbation theory or chiral perturbation theory 
which would yield a "mean field exponent" i' — 1/2. 
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independent of the interesting question what would be the true universahty 
calss of a second-order transition in two- flavor QCD, as seen at very large 
correlation length close to T^. 
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